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Abstract 

In this paper, we study the Nflation model, in which a collection of massive scalar fields drives 
the inflation simultaneously. We observe, when the number of fields is larger than the square of 
ratio of the Planck scale M p to the average value m of fields masses, the slow roll inflation region 
will disappear. This suggests that in order to make Nflation responsible for our observable universe, 
the number of fields driving the Nflation must be bounded by the above ratio. This result is also 
consistent with recent arguments from black hole physics. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



The multiple field inflation implemented by assisted inflation mechanism proposed by 



Liddle et.al 



l| relaxes many limits for the single field inflation models, and is being become 



a promising class of inflation models. There has been many studies on it 



Dimopoulos et.al J4] showed that the many axion fields predicted by string vacuum can 
be combined and lead to a radiatively stable inflation, called Nflation, which may be an 
interesting embedding of inflation in string theory. Then the detafled study was made by 
Easther and McAllister [5j for quite specific choices of initial conditions for the fields. In 
Nflation model, the spectral index of scalar perturbation is always redder than that of its 
corresponding single field, which is given numerically in and is showed analytically in 



a a. 



Recently, 



1 , the ratio of tensor to scalar has always same value as in the single field case 10|, the 



non-Gaussianity is quite small 
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12|. There was some further studies [1J 
In single field inflation model, the occurrence of inflation requires the value of field must 
be beyond the Planck scale, which can be obtained by imposing the slow roll condition upon 
the field. However, when the number of fields increases, this value will decrease rapidly, and 
can be far below the Planck scale, especially when the number of fields is quite large. This 
is a remarkable and interesting point of Nflation model. 

In addition, in single field inflation model, when the value of field increases up to some 
value, the quantum fluctuation of field will inevitably overwhelm its classical evolution along 
the potential. In this case, the inflaton field will undergo a kind of random walk, which 
will lead to the production of many new regions with different energy densities. In some 
regions, the field will wander down along its potential, so the classical variance dominates 
the evolution again and then inflation is able to cease when the field reaches its bottom. 
However, in another regions the field will fluctuate up and inflation will keep on endlessly. 

nn 

This so called slow roll eternal inflation [14], [15[ , has been studied by using the stochastic 



approach 

The critical value of field separating the field space into the slow roll inflation region and 
eternal inflation region can be obtained by requiring the change of classical rolling of field in 
unit of Hubble time equals to its quantum fluctuation. In the case of single field, this value 



1 see also different result for Nflation with small-field potential 
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is far larger than the Planck scale, and so the end value of slow roll inflation. However, in 
the Nflation model, it seems that when the number of fields is added, the total classical roll 
of fields is weaken, while its total quantum fluctuation is strengthened, which will lead to 
this critical value moves faster to some smaller one, which maybe bring a bounds for the 
number of fields participating in inflation. Here the 'value' for multiple field means what 
we take is the root of square sum of changes of all fields, because here all fields contribute 
inflation, and thus the trajectory is given by the radial motion in field space, Thus it is 
interesting to check this possibility. This will be done in this paper. In section II, we will 
study the case of Nflation with massive fields. Firstly we show a simply estimate for the 
bound of fields number by taking Nflation with equal mass fields as an example. Then we 
study a general case with mass distribution following Marcenko-Pastur law proposed by R. 
Easther and L. McAllister [5|], which further validates our result. In section III, we discuss 
the case of Nflation with 4 fields. The summary and discussion is given in the final. 



II. BOUND FOR N OF NFLATION 

In Nflation model, the fields are uncoupled and potential of each filed is Vi = |m|0f. 
The total change of all fields is determined by the radial motion in field space. In the slow 
roll approximation, we have 



A0 = /^(A^ 



H 



,2 



where A0, ~ ^ and <pi ~ ^ have been used, and the factor with order one has been 
neglected. In the meantime, the total quantum fluctuation of fields is 



'£«.) 2 -,/£(£> 2 



N 



( 2 ) 



where N is the number of fields, 5<pi ~ ^- has been used and the factor with order one has 
been neglected. By requiring A</> = S(j), we will obtain the critical point separating the slow 



roll inflation region and eternal inflation region, which is given by 

£™?#) - f E( m M 2 - ( 3 ) 

i / i 

In slow roll inflation region, the end of slow roll inflation requires ~ 1, which may be 
reduced to 

*?EW^E m ^) • ( 4 ) 

A. The case with equal masses 

Firstly, when the masses of all fields are equal, i.e. m, = m, and also for simplicity we 
take the values of all fields are also equal, i.e. fa = <fi. From Eqs. ([3]) and (jl]), we have 



1 M* 



Mp (6) 



<N 

respectively. Thus we see that the end point moves with -7=, which is slower than that of 
the critical point separating the slow roll inflation region and eternal inflation region, since 
the latter changes with ^74. This suggests that when we plot the lines of the end point and 
the critical point moving with respect to N, respectively, there must be a value for these two 
lines to cross. Beyond this value, the slow roll inflation region disappears as shown in Fig. 
CD This value can be obtained by taking both Eqs. ([5]) and ([6]) equal, which gives N ~ — §-. 
Thus to make Nflation responsible for our observable universe, the number iV of fields in 
Nflation model must satisfy 



M 



N<^ (7) 

since the existence of such a slow roll region is significant for solving the problems of standard 
cosmology and generating the primordial perturbation seeding large scale structures of our 
universe. It should be noted that in the case that the masses of all fields are equal, their 
field values are equal is not reality, because even if initially the values of all field are equal, 
they will also be unequal after several efolds due to the random walk of each field. However, 
this simplified analysis actually provides a simple estimate for the bound for N. In next 
subsection, we will validate this result in a general case. 
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FIG. 1: The figure of log change of the end point of slow roll inflation and the critical point 
separating the slow roll inflation region and eternal inflation region with respect to the number 
iV of fields in Nflation model with massive fields, in which m = 10~ 6 M p is taken. Three regions 
separated by both lines have been pointed out in figure. We can see that when N ~ 10 12 , the slow 
roll inflation region disappears. 

B. The case with mass distribution following Marcenko-Pastur law 

Then we will study a general case with mass distribution following Marcenko-Pastur law 
proposed by R. Easther and L. McAllister 5j, which appears for axions in string theory. The 
shape of the mass distribution of axions depends on the basic structure of the mass matrix, 
which is specified by the supergravity potential. In the simplest assumption, the mass 
matrix is a random matrix. When one diagonalize this matrix, the fields will be uncoupled 
with the mass spectrum given by the distribution of eigenvalues. This distribution of the 
eigenvalues can be characterized by Marcenko-Pastur law when the matrices are large. The 
mass distribution taken as the Marcenko-Pastur law is a function with respect to rh and 
(3, where m is the average value of the mass, i.e. < m 2 >= m 2 , and (3 is determined 
by the ratio of the number of axions to the dimension of the moduli space and a model 
dependent parameter, whose favourable value is expected to be about 0.5, see Refs. {5, 7\. 
In this case, the smallest and largest mass are given by m 2 = a = m 2 (l — \/j3) 2 and 
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m N = o = m 



1 + y/j3) 2 , respectively. 



In the slow roll approximation, the field value can be given by 

^-utoMt)}^, (8) 

where r(i) is the ratio of the value of the heaviest field at time t to its initial value to, 
r(t) = ^^p) • Then defining c = 21n ^" , the parts including (p 2 in the summation terms 
of Eqs. (jHJ) and (HI) can be replaced with 0|(t o ) exp [cmf]. When we ignore correlations 
between the mass distribution and the initial field distribution, we can straightly calculate 
their respective average values. By using power series expansions, the average value of 
exponential term exp [cm 2 ] can be written as 

rj 



< exp [cm 2 ] >= < > — . (9) 



The expectation value inside of summation in left hand side of Eq. (jHJ), can be expressed 
with Narayana numbers T(i,j) (see Eq. (6.14) in Ref. [5]) 

% 

<mf> = mV^TiiJ)?- 1 

3=1 

= m^F^l-it-i^P), (10) 

where 2 -Pi is hypergeometric function. Then Eq. ([9]) can be rewritten as 

00 j 
<exp[cm2] >=53m 2i 2 F 1 (l-i,-i,2,/3)^. (H) 

i=0 

Therefore the summations terms in Eq. ([3]) with the expectation values of initial conditions 
and of distribution of mass spectrum of fields are 

00 j 
J2 mlti = Nam 2 ^ m^F^-i, -i - 1, 2, $) C - (12) 

i i=0 

00 j 

]Tm^ 2 = iVam 4 ^mViH - 1, -i - 2, 2,/3)~ (13) 

i i=0 

where a =< <f>j(to) >. 

Thus the Eq. ((2]) with the help of Eqs. (fT2l and (1T31) becomes 

M 3 

a*—fcfi(t,0), (14) 
miv 2 
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FIG. 2: The figures of f\(t, f3) and /2 (£,/?) with respect to c with different (3. The average mass 



m = 10 M p is taken. 



where 



1/2 



Er=o^ 2 ^i(-j,-j-i,2,/?)^ 



■l 3/2 



(15) 



whose dependence on different c is plotted in Fig. [2j We see that fi(t,(3) is approximately 
a constant with order one for a wide range of c, i.e. different initial conditions and values of 
fields. Further, we can note that when all fields have equal values and masses, fi(t,f3) will 
have the value f±(t, 0) ~ 1 with a = 2 and rh = m. In this case Eq. (|14|) will be exactly 
same as Eq. ([5]). 

The Eq. (0} with the help of Eqs. (JT2]) and ([13]) becomes 



Ml 

-j*-Mt,l3), 



(16) 
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where 



l,-i-2,2,/ty 



2 ' 



(17) 



whose dependence on different c is also plotted in Fig. [2J We see that similar to fi(t,fi), 
/2 is also approximately a constant with order one for a wide range of c. When all 
fields have equal values and masses, fait, 0) — 1 with a = </> 2 . In this case Eq. ffT6]) will be 
exactly same as Eq. (J6]). 

Thus combining Eqs. (|14p and (fT6l) to cancel a, we have iV ~ ^f-, where the ratio of 
fi(t, (3) to f2(t, P) has been taken as roughly 1, which can be seen in Fig. El This is a point 
in which the slow roll inflation region will disappear. Thus to have a period of slow roll 
Nflation, the number of fields must be bounded by 



N < p 



rsj — o 

m 2 



(18) 



This result further validates the argument in previous subsection, only replace m with m. 



III. THE CASE OF NFLATION WITH cp 4 FIELDS 

It is interesting to further check whether there is similar bound for the field number of 
Nflation with 4 fields. Following the same steps as we did in the previous section, the 
critical point separating the slow roll inflation region and eternal inflation region and the 
end point of inflation can be given by 

-^FB A ^) 2 > ( 19 ) 
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M p 2 £(A^) 2 ^ (E A ^J > ( 2 °) 

respectively, where A, is the couple constant of the corresponding field and the factors with 
order one have been neglected. For simplicity, we take all = <fi and Aj = A, and thus have 



M p 



N 



(21) 
(22) 
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FIG. 3: The figure of log change of the end point of slow roll inflation and the critical point 
separating the slow roll inflation region and eternal inflation region with respect to the number N 
of fields in Nflation model with (ft 4 fields, in which A = 10~ 12 is taken. Three regions separated by 
both lines have been pointed out in figure. We can see that both lines are parallel, thus there is 
not the bound for N. 

Thus in this case the critical point and the end point approximately obey the same evolution 
with iV increased, which is plotted in Figj3j This suggests that there is not the bound for 
the number of fields imposed by the occurrence of slow roll inflation. 

This result seems unexpected. The reason leading to it may be that for 4 field, its 
effective mass is ~ X(ft 2 , which is changed with (ft, and its change in some sense sets off the 
fast moving of the critical point. When writing X(ft 2 = m 2 in Eq. (j2ip . one can find that the 
resulting equation will be the same as Eq. (jSJ). Thus combining it with Eq. (1221) . we will 
have the same result with Eq. (jTJ), which in turn suggests 



Thus for (ft A field, the bound relation between the mass m and N in Eq. (JTJ) for massive field 
is transferred to that between the field value (ft and N. The study with general case will be 
expected to have approximately the same result with Eq. ( 1231) . which is neglected here. 



X(ft' 




(23) 
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IV. SUMMARY AND DISCUSSION 



In this paper, we study the Nflation model, in which a collection of massive scalar fields 
drives simultaneously the inflation. We observe that with the increase of fields number, both 
the critical point separating the slow roll inflation region and eternal inflation region and the 
end point of slow roll inflation will move towards smaller average value of fields, however, at 
different rates. In general, the critical point moves faster than the end point, which leads to 
that the slow roll inflation region will be eat off by the eternal inflation region inch by inch. 
When the number of fields is enough large, i.e. N ~ M 2 /rh 2 , both points overlaps, which 
means that the slow roll inflation region completely disappears. In this sense, in order to 
make Nflation responsible for our observable universe, the field number driving the Nflation 
must be bounded by iV < M 2 /m 2 . 

Recently, it was shown that in theories with a large number N of fields with a mass scale 



20 



2l|, i.e. Nm 2 < M 2 which 



m, black hole physics imposes a bound between M p and TV 
is actually the same as the result showed here. This can be explained as follows. In general 
each field can contribute the factor ~ m 2 into the renormalization of the Planck mass, thus 
after the accidental cancellations are neglected, the net contribution of iV fields will be 
Nm 2 ~ M 2 . This indicates that with N massive fields there exists a gravitational cutoff, 

fl. When we focus on the 



beyond which the quantum gravity effect will become important 
inflation driven by N massive fields, we observe that the same cutoff will also appear in a 
similar sense, i.e. beyond this cutoff the quantum effect will be dominated, it is which that 
leads to the disappearance of the slow roll inflation region. Thus in this sense our result also 



justifies the bound of Ref. 



20] from a different point of view. In this case exactly the Plank 



mass should be replaced by renormalized one, i.e. M 2 , which includes the contributions 
of all massive fields for M 2 . However, it can be noted that M 2 is actually the same order 
as M 2 . Thus when the renormalization of M 2 is considered, our result is not qualitatively 
altered. 

However, this bound can not be applied to nearly massless scalar field. The reason is 
when the masses of fields are negligible, they will not appear in summation for fields in 
both sides of Eqs.Q and (jlj), which is actually also a reflection that the massless fields 
do not affect the motion of massive fields dominating the evolution of universe, while the 
perturbations summed in Eq.(2) are those along the trajectory of fields space, since the 
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massless fields only provide the entropy perturbations orthogonal to the trajectory, which 
thus are not considered. The same case can be also seen in argument of 20|, since the 
contribution of field to the renormalization of the Planck mass is proportional to m 2 , thus the 
net contribution leaded by all massless fields to the Planck mass may actually be neglected. 
Thus if there are some nearly massless fields and some massive fields with nearly same order, 
it should be that there is a bound iV < M 2 /m 2 , in which only massive fields are included 
in the definition of rh and N. In our example when the mass distribution is characterized 
by Marcenko-Pastur law, as has been used, in which the smallest and largest mass are given 
by m 2 = m 2 (l — \fj3) 2 and m 2 N = m 2 (l + \/j3) 2 , respectively, the masses of all fields are 
approximately in same order for (3 ~ 0.5. In this case it is natural that all fields need to be 
considered. 

The field A0 4 in essence is different from massive field m 2 (f) 2 . The former cor resp onds to 
have a running mass ~ A0 2 , which is dependent of <fi. In this case, following Ref. [20], the net 
contribution of iV fields to the renormalization of the Planck mass will be ~ NXcj) 2 ~ M 2 . 
Thus a same bound relation with Eq.( [23l) can be obtained, which is seemingly one between 
the field value and N. Eq.( }23]) can be written as 



< 



(24) 



r ~ VXN 

Note that for general A < 1, when there is a slow roll inflation region, Eq. f[2~4|) is always 
satisfied, since the inequality given by Eq. ffT9l is actually included in Eq. (l24l . Thus it seems 
that there is not a bound for N in the Nflation with A0 4 . 

In principle, the bound showed here seems be only valid for massive scalar fields. Further, 
whether there are similar bounds for other fields with various potentials remains open, and 
needs to be studied. Be that as it may, however, the result observed, that there may be 
a large N transition leaded by the quantum effect in inflation, may be interesting, which 
might have deep relations with other large N phenomena discussed, and thus is worth further 
explore. 
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